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Abstract 

The singularity for the big bang state can be represented using the 
generahzed anisotropic Friedmann equation, resulting in a system of 
differential equations in a central force field. We study the regulariz- 
ability of this singularity as a function of a parameter, the equation 
of state, w. We prove that for > 1 it is regularizable only for w sat- 
isfying relative prime number conditions, and for w < 1 it can always 
be regularized. This is done by using a McGehee transformation, usu- 
ally applied in the three and four-body problems. This transformation 
blows up the singularity into an invariant manifold. The relationship 
of this result to other cosmological models is briefly discussed. 



1 Introduction 

It is well known that the Kepler problem, which defines the motion of two 
bodies moving in a Newtonian gravitational inverse square central force field, 
can be regularized at the collision of the two bodies. Letting r represent the 
distance between the two bodies, then the regularization of collision implies 
that there exists a change of position and velocity coordinates, and time, 
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so that in the new coordinates and time the state corresponding to r = 
is well defined. Moreover, the flow of the transformed differential equations 
in position and velocity space is smooth(real analytic) in a neighborhood of 
collision, and also on the entire associated fixed energy surface. This implies 
that solutions can be smoothly extended across collision as a function of time. 
The Levi-Civita transformation is a classical example of a regularization that 
accomplishes this in the two-dimensional Kepler problem (Levi-Civita 1920). 
Physically, this means that the two bodies perform a smooth bounce at 
collision. Within the field of celestial mechanics, there are many other regu- 
larizations for binary collisions in the Kepler problem. These include the use 
of geodesic flow equivalent transformations, due to Moser and this author 
(Moser 1970; Belbruno 1977), and the Kustaanheimo-Stiefel transformation 
(Kustaanheimo and Stiefel 1965). 

Regularization transformations can be applied to problems other than the 
classical Newtonian inverse square gravitational force, to which they have 
been generally historically focused. The regularization of the collision of a 
zero mass particle, eg a photon, with a Schwarzschild black hole was de- 
scribed in a recent paper by F. Pretorius and this author (Belbruno and 
Pretorius 2011). This is accomplished by equivalently viewing the motion of 
the particle as being in an inverse forth power central force. The differential 
equations describing this motion can be transformed using a mapping due 
to R. McGehee (McGehee 1981), which blows up the collision state into an 
invariant manifold. This transformation describes the motion of the particle 
about the black hole is a way that eludicates the entire phase space in a sur- 
prisingly clear fashion. At the same time, the collision state of the particle 
with the black hole can be regularized. This turns out to be a 'branch regu- 
larization', where individual solutions have unique extensions across collision 
as a function of time, which, in position coordinates, are only continuous at 
collision and not differentiable there. The key point is that the extension 
across collision is uniquely deflned. This implies that when the particle col- 
lides with the black hole, it performs a continuous, but not differentiable, 
bounce to a unique extension, or bounce, from the black hole. 

In this paper we can gain interesting information about the dynamics 
in a neighborhood of the singularity due to the big bang using a McGehee 
transformation, following the same approach used in Belbruno and Preto- 
rius(2011). We flrst show that dynamical flow near the big bang singularity 
can be reduced to a central force fleld, when modeled by an anisotropic Fried- 
mann equation, under a number of assumptions. This force fleld is deflned 
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with inverse powers of Qj ; ^ 

> 5, where a is the expansion factor of the 
universe and a = corresponds to the big bang. We then apply the McGehee 
transformation to the central force field, yielding unique branch extensions 
of solutions through a = 0. 

We prove that the big bang can be branch regularized, provided that one 
of the parameters, the equation of state variable, w > 1, satisfies relative 
prime number conditions, which are a set of measure zero. In other words, 
for w > 1, except for a set of measure zero, the big bang is not branch 
regularizable. Under these assumptions, this implies that if we consider a 
universe prior to ours, for t > that collapsed to a = at time t = 0, 
then as a function of time, our universe for a(t) > 0, with t < 0, cannot 
be defined as a branch extension, from the previous universe for w > 1, 
except for a set of measure zero. However, for w < 1, the singularity can 
always be branch regularized. The equation of state, w, is a key parameter in 
cosmology, w = P/p, where P is the pressure and p is the energy density of 
the universe. The parameter w and relevance of if = 1 is explained further 
in Section HI It is noted that we have chosen the convention that t > 
prior to the big bang singularity and t < after the big bang. This is 
chosen for convenience to be consistent with the time convention used in 
McGehee(1981), since McGehee(1981) is referred to frequently. 

This result can be related to other cosmological models, which we do in 
Section |H It is important to note that our results are purely mathematical 
in nature and we are not making any assertions to the relevance of these 
results to reality. 

The paper is organized as follows: In Section |2]the problem is formulated. 
Section |3] gives the solution to the problem and the main results. Block 
regularization is discussed at the end of Section [31 In the final section, 
Section HJ we briefly discuss the relevance of our results to other cosmological 
models. 

2 Formulation of Problem 

To formulate the problem, we begin with the cosmological principle. This 
states that at any moment in time, the universe presents the same aspect 
from every point, except for local irregularities. This equivalently means that 
the universe is isotropic and homogeneous. 

Assume that we know the distribution of the total finite number of galax- 
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ies at time to of distance rj from a fixed origin O, i — 1,...,N where N 
is the total number of galaxies. Then the cosmological principle implies 
that ri{t) = a{t)ri{to), where t G 3? is time. a{t) is called the Friedmann- 
Robertson-Walker(FRW) scale factor. That is, the cosmological principle 
implies a uniform expansion or contraction. Hubble's law implies, 



where H = H{t) is the Hubble parameter, a = a(t), ■ = d/dt. H{to) is the 
Hubble constant, which is the current expansion rate of the universe. 

We consider the Priedmann equation that approximates cosmic evolution. 



where a = a{t) is normalized, so that for today(t = to < 0) a{to) = 1. p° is 
the present value of the energy density for component j, where m represents 
non-relativistic matter, r represents radiation and w represents an energy 
component with equation of state w. The last two terms on the right-hand 
side represent the spatial curvature and anisotropy, respectively. K is the 
curvature and a is the anisotropy parameter (see Garfinkle et al. 2008). 

We are interested in studying solutions where a — >■ 0. This is initially 
assumed to take place in a contracting phase of the universe, just prior to the 
big bang, where t > 0. This assumption is consistent with current alternative 
theories in cosmology to the standard big bang inflationary model, such as 
the cyclic or ekpyrotic models (Steinhardt and Turok 2002; Garfinkle et al. 
2008), where the universe undergoes a period of slow contraction with w > 1 
prior to the big bang. We are assuming that the big contraction (crunch) 
converges to a = 0. 

The goal of this paper is to understand under what conditions solutions 
a — a{t) can be continued, in a well defined manner, with i > 0, as a 
function of time to t = 0, which corresponds to the start of the big bang, 
and then for t < beyond the big bang. The way this will be done is to see 
under which conditions, the function a{t) can be regularized at t = 0. The 
type of regularization we are using is branch regularization, mentioned in the 
Introduction, where the trajectory for a{t) for t > is smooth(analytic) and 
can be uniquely extended to t < 0, which is continuous, at t = and analytic 
for ^ < (see McGehee 1981). Branch regularization yields a continuous well 
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defined unique bounce at t = for tlie given trajectory witli initial condition 
in tlie contracting pliase for t > 0. We explain how this can be accomplished 
in the next section, where branch regularization is carefully defined. 



3 Assumptions, Differential Equations, and 
Results 

Friedmann's equation, ([2]), gives rise to a second order differential equation 
for a{t) under some assumptions on the cosmology model we are using. 

Assumption A Let t = ti be the initial time for a{t) in the contracting 
phase of the universe prior to the big bang at t = 0, so that < ti ^ 1. 
Since it is assumed that a — >■ as t — )■ 0"*" (O"*" means t — )■ with t > 0), then 
a(ti) is assumed to be sufficiently small, say < a(ti) < 6 1. 

Assumption B w and K are constant, a, p^, p^, are positive constants. 

It is reasonable to fix these parameters in a neighborhood of a = 0. For 
example, in the ekyprotic model of the big bang, where w > 1, numerical 
computations show that these parameters generally converge to constant 
values (Garfinkle et al. 2008). 

Under Assumption B, we differentiate ([T]) with respect to t. We first 
substitute into([I]), which yields 

3 V a a a^(i+w)-2 J ^ ^ 

Applying d/dt to both sides of ([3]) yields 

We can divide through both sides of (jl]) by d which is nonzero. This follows 
from ([3]) which implies that as a — )■ 0, |d| — t- oo. After rearranging terms, we 
obtain, 

•■ = _ , M , P^[3(l + ^)-2] 
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The differential equation given by ioi a = a{t) plays a key role in our 
analysis. It defines the behavior of a{t) as a function of w. We will consider 
three cases: w < 1, w = 1, and w > 1. 

The term 'branch regularization', mentioned in the Introduction, is for- 
mally defined as in McGehee(1981). We define it for a general first order 
system of differential equations, 

X = F(x) (6) 

where • = (i/cit, x G -R", n > 2. Assume F is a real analytic vector field on 
an open set U <Z BJ^. 

We first define the term 'singularity'. 

Definition 1 Let x(t) be a solution for dH]), defined in f/ for 6 < t < a, with 
initial condition x(to), —oo < b < to < a < oo. {b,a) is assumed to be the 
maximal interval the solution can be extended. If — oo < b then the solution 
is said to begin at t = b and if a < oo, then the solution is said to end at 
t = a. In either case, t* = a or t* = b is said to be a singularity for x(i(:). 

The terms 'branch extension' and 'branch regularization' are defined in 
Definitions 2,3, respectively. 

Definition 2 Let Xi(t),X2(t) be two solutions of IQ. Suppose that Xi 
ends in a singularity at t* and X2 begins in the singularity at the same time. 
Assume there is a multivalued analytic function having a branch at t*, and 
extending both Xi,X2, where we regard t as complex. Then Xi is said to 
be a branch extension of X2 and X2 is said to be branch extension of Xi. 

Definition 3 A solution X of ([HD which either begins or ends in a singu- 
larity at time t* is called branch regularizable at t* if it has a unique branch 
extension at t*. 

It is important to note that in Definitions 2,3, time t is assumed to be 
complex. The solutions Xi,X2 are real analytic, taking on real values for 
real t. However, they are extensions of each other through complex t. This 
definiton of branch regularization is well defined for one degree of freedom, 
n=l, as is seen, using time complex variable. 

We are interested in obtaining branch extensions for ([5]) at t = t* = 0. 
([5]) is written as a first order system analogous to ([6]) by setting x = (a, a). 
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We need only find a branch extension for for a{t) since that automatically 
yields one for a(t). 

Branch regularization is defined along trajectories of solutions. It's main 
feature is that it yields unique extensions at singular points whose solu- 
tions are not defined in the original coordinates. Branch regularization at 
a = for ([5]) will be proven by a change of coordinates and time, a, a, t, 
respectively, with an explicit transformation, to new coordinates, r,v,s, re- 
spectively, where the state a = 0,\a\ = oo is well defined, corresponding to 
r = 0,v = ±y/2. In this way, the singular state at a = is resolved. The 
solution, a{t), is then well defined and continuous at t = 0. This is obtained 
through an explicit formula. 

It is remarked that although McGehee's work (McGehee 1981) is formulated 
in general for two degrees of freedom, it can be reduced to one degree of 
freedom, as we have applied in this paper, in a straight forward manner, 
when extending a solution a{t) from t > to t < 0. The proof of a unique 
extension in McGehee(1981), in general, requires knowing the explicit de- 
pendence of the solution, in this case a{t), as a function of t near t = by 
transforming a and d to regularized coordinates, r, v, using the original time 
coordinate t. In the regularized coordinates, the explicit form of the solution 
r(t) can be obtained using the stable manifold theorem, where uniqueness is 
established using both r, v. This yields a unique branch in these coordinates. 
Transforming from r to a yields a unique branch for a{t). This is described 
in more detail in the proof of Theorem 5. 

We are now able to state our results. 

Theorem 1 

(a.) If |w| < 1, then every solution a — )■ can be branch regularized at 
a = 0. 

(b.) If w > 1, then a branch regularization of a = can be done on a 
countable set (measure zero), Q^, of w values, or equivalently, a branch reg- 
ularization of a = does not exist for almost all values of w (i.e. outside the 
set Qw)- The set is given by 

Q^ = {w = ^-1 \{p,q) e p}, 
where p, q are positive integers and 
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P = {{p^ < P < (Pj q) relatively prime, q odd }. 

The condition of regularizability with w G Qw is both necessary and suffi- 
cient. 

We now show how this theorem is proven in a series of theorems. 

A preliminary constant time scaling is done from t to r (this use of r is not 
to be confused with its earlier use in the symbol, p°), 

r2(2)-i/V-ir, if < 1, 

t = I 2C-^/^T, a w = l, (7) 

[{{A/3)7rGpl)-'/\ if ^>1, 

where 

C = 2a^ + (16/3)7rGp° . 

By Assumption B, this time scaling is well defined, where C > 0. 
As a — )■ 0, we can take a sufficiently small so that the term with the largest 
power of in the right hand side of (|5]) dominates the behavior of d. In 
fact, the following theorem is directly obtained. 

Theorem 2 In the scaled time variable, r, given by ([7]), and where ' = 
d/dr, the differential equation (jS]), can be reduced to the following: 

If w < 1, 

a" = -l-/i(a), (8) 

where 

/^(a) = aia"^ + a2a~^ + asa"^^"''!^ (9) 

K,i{w) = 3{l — w) > is a constant, are constants, and Ki — t- 0+ as w t 1- If 
we set c = (4/3)7rG', then ai = 2a^'^cp^, a2 = Aa~'^cp'^, = 2cr"^c(4 — fi;i)p° . 
Thus the values of the powers of of /i are less than that of the ffist term 
on the right side of ([8]). 

If w = 1 

a" = -4 -/2(a), (10) 
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where 

/2(a) = M~' + M-', (11) 

and bk are constants; 61 = 4cC~^p[^, 62 = 8cC~^p°. 
If w > 1, 

-" = -^-/3(«X (12) 

where 

/3(a) = cia"^ + C2a^^ + caa"^, (13) 

where K2 = 3(w — 1) > is a constant, are constants, and K2 — ?■ 0"*" as 
w I 1; ci = pZiPw~^-,C2 = 2p°p° ~\ C3 = 2cr^p° ~\ Thus, the values of the 
powers of of /a are less than that of the first term on the right side of 

The coefficients ak,bk, Ck are well defined by Assumption B. 



The differential equations, ([H]), ( |TOl) ( |T2i) . can be written as 

where \^ is a central force potential function, 

V = ^ + fia), (15) 



where 



a 



4, if w < 1, 



3(1 + u;) -2, if w > 1. 
f{a) = 0{a~~y),y < a is a sum of terms of the form a"^'', pk < a. 
More exactly, 

{aia^i + \a2a-'^ + ^^-(4-^1)^ if w < 1, 
hia~^ + |62a"^, if w = 1, 

cia"^ + |c2a"^ + \cza~^, if w > 1, 

where it is noted that the term 03/(4 — /ti) is well defined for ki 
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The new variables are now a,P = a', t and the new differential equation is 
(IT^ . This can equivalently be written as a first order Hamiltonian system 
of one degree of freedom, 




P' 



dH 
da 



(16) 



where, 



H = -P^- Via 



(17) 
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H is an integral of the motion, so that each trajectory X{t) = {a{r), P{t)) 
lies on the energy manifold 



for some constant h G R . S is a one-dimensional curve, h is uniquely 
determined from the initial value of X{t). It is seen that on each manifold, 
S, |P| — 7- oo as a — 0. 

It is noted that by the time scaling ([7]), multiplying t by a positive constant, 
A, T = At, the results we obtain wrt r are immediately obtained for t, where 
r — i- is equivalent to t — )■ 0. Thus, we will present the results from this 
point on for (IT^ without loss of generality. 

Set 7 = (1 + /3)~^ and /3 = a/2, We apply a transformation that smooths the 
singularity at a = so that |P| = oo gets mapped into a well defined set. 

Theorem 3 The transformation Tm : a, P ^ r,v, of fll4p . or equivalently 
(ITB]) . which smooths the singularity at a = 0, is given by 



{a,P\H{a,P) = h}. 



(18) 



a = r 



P 



r 



(19) 



where, if w < 1, 



a = 4, 



/3 = 2, 



7 = 1/3, 



(20) 



and if w > 1, 




This case implies that < 7 < 1/3 and /3 > 2. 
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The Hamiltonian energy (fT7|) gets mapped into the transformed energy equa- 
tion, 

t;2-2 = 2/ir°^ + 2r"^^(r), (22) 

where 

g(r) = f{r^). 



The proof of this is obtained by noting that the transformation Tm is di- 
rectly obtained from a more general map for two-degrees of freedom in McGe- 
hee(1981)(see Equation 4.1, where w = 0,9 = 0, and Equation 4.7). It is sim- 
ilarly verified that ( 122|) is obtained(see Equation 4.3 in McGehee(1981)(see 
also Belbruno and Pretorius 2011)). In our case, we have the additional 
term /(a) which gives rise to the term 2r°''^g(r) = 0{r^"~y^"') in ( 122|) . where 
a — y > 0, and 7 > 0. More exactly, 

{air + ia2r2/3 + -^r'^i/s^ jf w < 1, 

hr + i62r2/^ ' if = 1, (23) 

C^^7(3+.C2) ^ i^^^7(2+«2) ^ lc3rT"'2, if W>1. 

Thus, the energy manifold S is mapped into 

M = {r,v\v^ ~2 = 2hr'''' + G{r)} (24) 

Since a = is mapped into r = 0, and G'(O) = 0, then this defines the 
'collision manifold' N (see McGehee 1981), which represents the big bang 
singularity, 

N = {(r,f) eM\r = 0,v = ±V2}. 

Thus, |P| = 00 is mapped into v = ±\/2. A set analogous to this is defined in 
Belbruno and Pretorius(2011) for the problem of motion about a black hole, 
using the two-degree of freedom version of Tm- N is an invariant manifold 
for the transformed fiow in a new time variable we define below. 

The transformation Tm smooths the singularity at a = in the sense that 
\P\ = 00 is mapped into a finite set of points. To see if a solution a = a(r) can 
be extended through collision as a function of r, it is necessary to examine 
the transformed differential equations. 



11 



Set 

Qy = {i = ^ \ip,q)^p-} (25) 

Adapting a theorem in McGehee(1981) of two degrees of freedom to one de- 
gree of freedom, and in particular to our system f|T6l) . and also using Equation 
( l28l) . proven to exist below in the proof of Theorem 5, we have 

Theorem 4 (McGehee) Every singular solution a — > to (fT6l) can be 
branch regularized if and only if 7 G Q^. 

This implies that if 7 ^ (^7, then no singular solution can be branch regular- 
ized, and if 7 G Q-y, then every singular solution can be branch regularized. 

The relative prime number condition in Theorem 4 is a key condition for this 
analysis. It results from a more general theorem in McGehee(1981), proven 
for the more general case of two-degrees of freedom for ( fT6l) . where, more 
generally, a = (ai,ai), P = (^1,^2), and in ([IT]), = |Pp, a" = |a|". 
However, it can be applied to our case of a single degree of freedom in a 
straight forward manner, as is verified. 

Before we look at the transformed flow and the explicit form of the solution 
a(r), we discuss the set Q^,. 

The prime number condition that 7 G Q.y is equivalent to a prime number 
condition on a and w. This follows from the definition of 7 = (1 + (a/2))~^ 
which implies, 

a = 2(7-^ - 1) 

This yields the condition that 

a G Qa, 

where 

Q^ = {a = 2ir'-l)\ieQ^}. 

We also obtain a prime number condition on w. This follows from the con- 
dition q; = 3(1+w) — 2in Theorem 3 for w > 1, which implies, 

a — 1 
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This implies that for w > 1, using the prime number condition in Qa for a, 

2 

w e Qi, = {w = -7"^ - l|7 G Qy}. 

that was given in Theorem 1. By Theorem 4, branch regularizabihty is 
guaranteed. 

It is noted that for the case of w < 1 in Theorem 3, j = 1/3 E Q-y guaran- 
teeing branch regularizabihty. 

This proves Theorem 1. 

Why is the condition that 7 G (^7 required? 

Let a = Fi{t) be the singular solution branch fi{t) in the scaled time 
variable r, where Fi — )■ as r — )■ 0^, ( t = is mapped into r = 0). We 
prove below, using the stable manifold theorem, what the exact form of Fi 
is in terms of powers of r which is given by fl28|) . This yields an expression 
for Fi{t) where the term t'^ appears. If 7 were irrational or of the form 
p/q, p > 0,q > with q even, p,q relativly prime, then r'^ would not be 
defined for r < 0, and Fi{t), r > 0, could not be extended to a branch F2{t) 
for r < 0. This gives an obstruction to the continuation of a solution a(r) 
through r = due to the appearance of imaginary time. 

We now state the result on the transformed flow and the form of the solution 
a(r). 



Theorem 5 Tm, together with the transformation of time from r to s, 
dr = rds, maps f|T6|) into the first order system 



nv n'i) 

- = (/3 + l)rt;, - = p(v'-2)-Gir), (26) 



defined on M, where 



2 '^1 

air + a2rs + a^r^ , if w < 1, 



G{r) = { bir + b2r-^, if w = 1, (27) 

C^^«27+3 _^ c2r'^27+2 _^ C3r^21^ if w >1, 
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The solution for a = Fi(r) — ?■ as r — t- 0^ can be written in explicit form as 

a = Fi(r) = rT^(r'^\r'^2), r > 0, (28) 

{ui,U2) = (07/4, Ki/3), (1/3, 1), (fi;27,7) for w < 1, w; = 1, w > 1, respec- 
tively, where \1/(X, Y) is a uniquely determined real analytic function of X, Y 
for (X, Y) in a sufficiently small neighborhood of (0, 0), 

vI/(0,0) = (v^(/3 + l)r. 

If 7 G Q.y, then the unique branch extension of Fi{t) to r < 0, where |r| <^ 1, 
is given by 

a = F2{t) = Fi{r) > if p is even, 
a = F2{t) = -Fi(r) > if p is odd, 
where Fi(0) = ^2(0) = ^(0,0). 

Proof of Theorem 5 

The explicit solution a = a(r), given by fl28l) . is obtained by considering f l22|) 
and (|23|) . in the original time variable r. This implies, 

^ = (1 + /3)(2 + 2hr''^ + G(r))i/2. (29) 
ar 

It is noted that a = 4, 7 = 1/3 for ti? < 1 and a = (4 + /t2), 7 = (1 + a;/2)~-^ 
for w > 1, fi;2 > 0, < 7 < 1/3. It is verified that for w < 1, the right 
hand side of of ( 129|) can be written as a analytic function fi{X,Y), where 
X = r'^^I^^Y = r*^^/^; for w = 1, it can be written as an analytic function 
f2{X,Y), where X = r^/^^Y = r; and for w > 1, it can be written as an 
analytic function f2{X,Y), X = r'^'^'^,Y = r^ . The functions /^(X, F) are 
real analytic for (X, y) sufficiently near (0,0) (to insure the square root is 
well defined) where A (0,0) = ^2(1 + /3). More exactly, 

/i = (1 + /3)(2 + 2hX^ + aiX^ + (1/2)02X2 + 03(^1 - A)-^Yf'^, 
/2 = (1 + /3)(2 + 2hX^ + hiY + h2X^f'\ 

/3 = (1 + I3){2 + 2/iXy^ + ciXY^ + (l/2)c2Xr2 + (l/4)X)i/l 
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These functions can all be expanded as a convergent Taylor series about 
(X,Y) = (0,0) with a positive radius of convergence. Thus, (12^ can be 
written as 

dv 

_ = ^(r-\r-^), (30) 

where uji{w) > 0, 102(10) > 0, and ^-"(0,0) = -\/2(l + (3), J^(X, F) is real ana- 
lytic for (X, Y) sufficiently near (0, 0). ((wi, W2) = («7/4, ki/3), (1/3, 1), {^2'^, 7) 
for w<l,w = l,w>l, respectively.) We can now apply the stable manifold 
theorem to ( 130|) in an analogous manner that was done in McGehee(1981)(see 
Lemma 5.7). This implies that there exists an analytic solution, 

r{T) = Tg{T^\T^^), (31) 

for r in a sufficiently small neighborhood of 0, Q{0, 0) = a/2(1 + /3). Setting 
a = r'^, we obtain 

and setting \1/ = yields the solution in Theorem 5. It is necessary to prove 
that Ui, U2 insure that t'^'=, = 1, 2, is not imaginary for r < 0. This follows 
immediately for w = 1, since Ui = l/3,Ci;2 = 1- For w < 1, Ui = 07/4 = 7 
and UJ2 = = l-w. But w = (2/3)(p/g)-l so that = 2[l-(l/3)(p/g)]. 
Finally, for w > 1, ijJ2 = ■y and Ui = K27 = 3(ty — 1)7 = 3(2[(l/3)(p/g) — 
l])7 = 6[(l/3)(p/g)-l]7. 

The system (I26p of differential equations is obtained by transforming 0161) 
using f|T9|) . together with the transformation of time from r to s, which is 
verified after some algebraic simplification. 

The concludes the proof of Theorem 5. 

Remark It is noted that the key result in Theorem 5 is Equation (1281) for 
r > as r — )■ O"*" and its continuation, F2{t), for r < as r — )■ 0~, in the 
original physical coordinates a, r. The flow in the regularized coordinates 
r,v,s given by Equation fl26|) . described below, is not required to conclude 
the extension of the solution. However, fl26|) is necessary to understand the 
nature of the flow near a = 0. 

Thus, by Theorem 5, a(r) = -Fi(r) — t- as r — )■ 0+, and at r = 0, 
Fi(0) = 0, then as r passes to negative values of r, a(r) = F2{t) increases. 
This motion of a provides a universe bounce through a = through a unique 
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branch extension. This motion as a function of r can be done as function of 
t through the constant scahng ([7]). By Theorem 4, the condition 7 G (^7 is 
required. The bounce is nondifferentiable at r = 0, or t = 0, which is seen by 
differentiation of Fi with respect to r. This introduces the term t'^~^ that 
multiphes Fi. Since 7 < 1, then r = is undefined. 

Geometry of flow in branch regularized coordinates 

Lemma 1 is an invariant manifold for the flow of (l26l) . As a trajectory 
approaches as r — )■ 0"*", then this implies s — )■ 00. 

This is proven from the definition of A^, which implies r = and v = ±\/2. 
These are rest points for ( |26l) . Thus, A^ is an invariant manifold set for (l26ll . 
This implies that in the time variable s, s — t- 00 as r — )■ 0^. 

It is noted that the vector field defined by the system of differential equa- 
tions fl26|) is not differentiable at r = 0, with respect to r, due to the term 
G{r). This is seen by noting that dO/dr is singular at r = 0. It is, how- 
ever, locally Lipschitz continuous in a neighborhood of r = 0, guaranteeing 
uniqueness. It is real analytic with respect to v aX v = ±v^ and r = 0. 

It is remarked that this branch regularization can be viewed as a mini- 
mal type regularization where the flow in the regularized coordinates can be 
smoothed to the extent where the collision set results in rest points. This 
yields unique branch extensions in physical space which are continuous but 
not differentiable. In other problems, eg the Newtonian inverse square cen- 
tral force problem, collision can be regularized so that in the regularized 
coordinates, the flow can be extended through the transformed collision set 
in finite time, where the flow is smooth in a full neighborhood of collision. 
These are 'global' regularizations, where, moreover, the flow is smooth on the 
entire energy manifold. Examples of these are mentioned in the Introduction, 
eg the Levi-Civita regularization or the geodesic equivalent flow maps (Levi- 
Civita 1920; Moser 1970; Belbruno 1977) (see also Belbruno 2004). These 
stronger regularizations also imply that in the original coordinates for the 
Newtonian inverse square central force law, each collision solution can be ex- 
tended through collision in a bounce as in the weaker branch regularization. 
The difference is that the global regularization yields a smoother behavior of 
the flow near collision than in the branch regularization. 

It is important to mention that another type of regularization of interest 
is 'block regularization'. In this case the flow of the differential equations 
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is viewed in a compact neighborhood of the singular set, called an isolating 
block. In one version of this, due to C. Conley and R. Easton (Conley 
and Easton 1971), the singular set is an invariant set where the vector field 
is defined. In another variation of this, due to R. Easton (Easton 1971) 
the vector field need not be defined on the singular set. In either case, 
the flow is studied near the singular set within the isolating block. The 
flow is called block regularizable, in a deleted neighborhood of the singular 
set, if it can be shown to be diffeomorphic to the trivial parallel flow. It 
is proven in McGehee(1981), in the case of a pure power law, a~", given 
by (ITSl) with f{a) = 0, that the flow is block regularizable if and only if 
P = 1 — n^^,n = 1,2, 3,.... This is done by showing an equivalence with the 
flow resulting from applying a Levi-Civita type regularizing transformation. 
A similar result may true in our case when /(a) is included. However, in this 
paper, (3 > 2. This result seems to suggest that when /(a) is included, the 
flow may not be block regularizable, although this is not proven. 

Block regularization does not require looking at the flow at the singularity, 
but near it. If it is block regularizable, then this implies a degree of uniformity 
of the flow near the singularity. One could make the rough analogy that the 
flow could be combed. Physically, this would imply that the dynamics of the 
system near the big bang, for example, would be fairly well behaved. If it 
were not block regularizable, then the behavior of the flow near the big bang 
would not be uniform in the sense that it would not be equivalent to a trivial 
parallel flow. 

It is noted that branch regularization does not imply block regularization, 
and although block regularization does imply branch regularization in the 
case of a pure power law, as is proven in McGehee(1981), it is not necessarily 
the case in general that block regularization implies branch regularization. 



4 Relevance of Results 

The big bang singularity is not understood in reality. The big bang is thought 
to have occurred approximately 13.75 billion years ago. The initial rapid ex- 
pansion of the universe after this state, modeled by the theory of inflation, is 
estimated to have occurred approximately within 10~^^ and 10~^^ seconds, 
with an expansion factor of minimally 10''® in volume. The initial tempera- 
ture near this period is estimated to have been 10^^ degrees Kelvin, beyond 
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our ability to experimentally reach in particle accelerators. There are a 
number of cosmological theories that try to understand the big bang. In this 
paper we have taken a purely mathematical approach, using a Friedmann 
model which is generally used to understand the expansion and evolution 
of the universe, and not the big bang singularity itself. We have applied 
a McGehee regularization transformation to this model, to try and under- 
stand the dynamics in a neighborhood of the big bang singularity, given by 
a = 0. This mathematical analysis does provide information on the flow of 
the differential equations we are using in this neighborhood, derived from 
the Friedmann model, yielding conditions on which solutions can be mathe- 
matically extended through the singularity. The existence of such solutions 
is not intended to be applicable to reality. Our analysis does not attempt to 
explain any of the physics or nature of this state. 

It is instructive to make a very brief comparison of our methodology 
and results to other cosmological models. For a comparison, we will use 
the Friedmann equation given by ([2]), which results from general relativity. 
We will also use the equation of state w, since that plays a key role in our 
results, and also for the sake of simplicity. The parameter w was defined 
in the introduction. For a static universe, w = —1/3. When the universe 
expands, w < —1/3. In Theorem 1, w = 1 is a transition case between 
extensions of solutions through a = that can be done, in general, for w < 1 
and for w > 1 where it is necessary that w G Qw As is described in Erickson 
et al.(2004), for w < 1, the anisotropy and curvature terms in the Einstein 
field equations grow rapidly and become dominant as a — )■ 0, that can cause 
chaotic behavior that may destabilize the contraction process. However, if 
w > 1, these terms remain negligible as compared to the energy density, p. 
In this sense, w = 1 represents a transition. This is interesting to compare to 
the role w = 1 plays as being a transition for conditions required for branch 
regularization. 

The cosmological theories we will be comparing too model the evolution 
of the universe using a Freidmann equation, however, their modeling of the 
big bang singularity are all quite different. It is important to note that 
there are many other factors that can be used in making such comparisons 
between theories, however, since we are not attempting to understand the 
big bang, the Friedmann equation for the evolutionary behavior plays a key 
role. Further analysis beyond this comparison is beyond the scope of this 
paper. There are many cosmology theories, and we are considering a subset 
for comparison. 
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In Gott and Li(1998), J. R. Gott and L. X. Lin discuss the possibility that 
our universe could expand from a = 0, then stop expanding and eventually 
contract into a 'big crunch' again to a = 0. This could repeat cyclically. 
In the contraction phase immediately immediately before a = 0, if w > 1, 
then by our Friedmann model, Theorem 1 would imply that for an extension 
through a = 0, it is necessary that w G Qw, a set of measure zero. However, 
if, as a — )■ 0, w < 1, then branch extension can be done in general. However, 
in this case, chaotic behavior can occur. 

P. Steinhardt, N. Turok, have put forward the ekpyrotic theory of the 
big bang (Steinhardt et al. 2001). It is a cyclic universe theory that uses 
string theory. In this theory, the big bang is hypothesized to have occured 
cyclically by the collision of branes, which are subspaces derived from string 
theory. The resulting expansion and contraction is modeled using a Fried- 
mann equation that is given by ([2]). In this theory, w > 1 and — )■ 1 as 
a — )■ (see Lehners 2009; Steinhardt and Turok 2005). This yields an inter- 
esting interpretation of results when compared to our methodology: If the 
convergence of if to 1 is within the tolerances required by physical consid- 
erations, including quantum fluctuations, yielding w = 1 in a neighborhood 
of a = 0, which is claimed, then an extension through a = is automatic 
by Theorem 1. However, if w > 1 is maintained for a > in the limit, no 
matter how small, then w would need to satisfy relative prime number val- 
ues on the set prior to a = to be consistent with Theorem 1, placing a 
restriction that w G Qw This condition, however, may be meaningless, since 
when a goes below the Planck scale, the laws of physics are no longer well 
defined. It is important to note that although there are issues for the ekpy- 
rotic/cyclic models for w > 1 m our analysis, quantum corrections, for any 
value of still allow an effective leading order description of the evolution 
of the universe in terms of the scale factor, a, governed by the Friedmann 
equation. 

There are other cosmological theories that use string theory. For exam- 
ple, another is string gas cosmology (Brandenberger 2004) put forth by R. 
Brandenberger and C. Vafa. This is a cosmological theory that incorporates 
a gas of strings within the four-dimensional space-time of the early stages 
of the inflation of the universe. The additional dimensions given by strings, 
offers a way to explain aspects of the universe in a way that is different than 
just looking at a pure inflationary model. The reduction of the evolution 
of the universe based on a string gas theory to a Friedmann equation was 
carried out by A. Kamenshchik and I. Khalatnikov (see Kamenshchik and 
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Khalatnikov 2011) by viewing a string g perfect fluid with w = —1/3. 

The Friedmann equation they obtain is different from ([2]), so that we would 
have to repeat our analysis with this equation to make conclusions on branch 
regularizability. Other theories using string theory include, brane cosmology 
and string landscape, which are beyond the scope of this paper. 

Another theory of the big bang we consider is by J. Hartle and S. Hawk- 
ing (Hartle and Hawking 1983). Their theory models the big bang using a 
quantum mechanical framework, where our universe is viewed as a particle 
and the big bang is modeled as a wave function. The evolution of the uni- 
verse after the big bang uses general relativity. The Hartle-Hawking theory 
assumes that time does not exist prior to the big bang singularity. This con- 
dition is inconsistent with the methodology of our results since we require 
a continuous time flow through the big bang. However, another theory of 
interest of the big bang using quantum mechanical ideas and quantized grav- 
ity, is called loop quantum cosmology (Thiemann 2007). It is based on loop 
quantum gravity, which is a quantization of general relativity, developed by 
A. Ashtekar. This theory gives rise to a cyclical universe theory where time 
can be continued though the big bang. This approach can be reduced to a 
modified Friedmann equation (Singh 2005). It is shown in Singh(2005) that 
w G [—1,1]. Within our methodology, these values of w this would auto- 
matically allow an extension of solutions through the big bang singularity. 
However, the Friedmann equation used in this theory is a modification of (12]), 
and we would have to repeat our analysis for this equation to make definitive 
conclusions. 

It is noted that Theorem 1 implies a 'unique' branch extension from a pre- 
vious universe to ours. As noted in the previous section, uniqueness results 
from Libschitz continuity of fl2^ . A modification of the original Friedmann 
equation could result in nonuniqueness of solutions at a = for the trans- 
formed differential equations. Also, we have fixed a number of variables as 
a — )• 0, such as w, and the curvature K, among others. Although fixing 
these parameters seems consistent with numerical simulations in Garfinkle 
et al.(2008), it would be interesting to study the approach taken in this paper 
if they were allowed to vary. This is a topic for further study. 
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